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We present charged particle densities as a function of pseudorapidity and collision 
centrality for the 197 Au+ 197 Au reaction at y / s/v r /v=130 GeV. An integral charged 
particle multiplicity of 3860±300 is found for the 5% most central events within 
the pseudorapidity range —4.7 < r/ < 4.7. At mid-rapidity an enhancement in the 
particle yields per participant nucleon pair is observed for central events. Near to 
the beam rapidity, a scaling of the particle yields consistent with the "limiting 
fragmentation" picture is observed. Our results are compared to other recent exper- 
imental and theoretical discussions of charged particle densities in ultra-relativistic 
heavy-ion collisions. 
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Multiplicity distributions of emitted charged particles provide a fundamen- 
tal measure of the ultra-relativistic collisions now accessible experimentally 
using the Relativistic Heavy- Ion Collider (RHIC). The particle densities are 
sensitive to the relative contribution of "soft" processes, involving the longer 
length scales associated with non-perturbative QCD mechanisms, and "hard" , 
partonic processes [1,2]. The total number of charged particles and the an- 
gular dependence of the charged particle distribution is expected to depend 
markedly on the amount of hadronic rescattering, the degree of chemical and 
thermal equilibration, and the role of subnucleonic processes. 

In this Letter we present charged particle distributions in dN c h/drj with —4.7 < 
7] < 4.7 for 197 Au + 197 Au collisions at v /^7=130 GeV. The pseudora- 
pidity variable rj can be expressed in terms of the scattering angle 6, with 
7] = — ln[tan(#/2)]. The data were measured at the RHIC facility using com- 
plementary subsystems of the BRAHMS experiment. We have also investi- 
gated the systematics of the pseudorapidity distributions as a function of 
collision centrality. The collision centrality can be related to the number of 
participant nucleons in the reaction, thus allowing for comparison of dif- 
ferent nuclear systems based on simple nucleon-nucleon superposition mod- 
els [3], as has been done, for example, at lower SPS energies for the Pb+Pb 
system [4,5]. Our results complement those recently reported by the PHO- 
BOS [6,7] and PHENIX [8] collaborations by establishing global dN ch /dr] sys- 
tematics through a very different experimental arrangement than used in those 
experiments. 

The BRAHMS detector system [9] consists of forward- and mid-rapidity mag- 
netic spectrometers, which allow for the determination of charged particle 
properties over a wide rapidity and momentum range, and a number of global 
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detectors employed to characterize the general features of the reaction, such as 
the overall charged particle multiplicity and the flux of spectator neutrons at 
small angles. The current analysis is based primarily on the global detectors 
(see Fig. 1 insert), including the Multiplicity Array (MA), the Beam-Beam 
counter arrays (BBC), and the Zero-Degree Calorimeters (ZDC) [10]. Also, 
the front time-projection chamber (TPM1) of the mid-rapidity spectrometer 
arm, which is an integral part of this spectrometer, provides independent con- 
firmation of the particle densities at selected angles. 

The MA is used to determine charged particle densities in the pseudorapidity 
range of —2.2 < rj < 2.2 for collisions at the nominal vertex. This pseudora- 
pidity coverage is increased to —3.0 < rj < 3.0 by using an extended range 
of vertex locations. For the present data, the collision vertex distribution was 
found to be roughly Gaussian in shape with a=70 cm. The MA consists of 
two layers of detectors arranged as six-sided coaxial barrels about the beam 
axis. The inner barrel (SiMA) is comprised of 25, 4 cmx6 cmx300 fim Si 
strip detectors located 5.3 cm from the nominal beam axis, with each detector 
segmented into 7 strips with a 0.86 cm pitch. For the SiMA, three sides were 
populated with six wafers each, covering 42 segments in pseudorapidity, one 
side was populated with five wafers, and the last two sides were populated with 
one wafer, each. The outer barrel (TMA) consists of 38, 12 cmxl2 cmx0.5 cm 
plastic scintillator "tile" detectors [11] located 13.9 cm from the beam axis. 
The scintillator tile array had four rows fully populated with 8 detectors, each, 
one row had 4 detectors, and the last row had 2 detectors. 

The single-particle response of the SiMA and TMA elements was calibrated 
by selecting peripheral collision events where energy peaks corresponding to 
the passage of single particles, mostly pions, through the detectors could be 
compared to the calculated energies for these peaks based on GEANT simula- 
tions [12] and assuming nominal detector thicknesses. The uncertainty in the 
single-particle response is estimated as ^5% based on the scatter of measure- 
ments for a given detector element for different vertex locations and based on 
the uncertainty in the GEANT simulations. 

For more central collisions at RHIC, the modest segmentation of the MA re- 
sults in multiple particles passing through individual detector elements. The 
number of particles passing through each element is determined by dividing 
the total energy observed in that element by the corresponding average energy 
loss for a single particle, primarily either a pion, kaon, or proton, as determined 
by GEANT simulations. Near to mid-rapidity and for central collisions an av- 
erage of 10 particles pass through the Si elements and 61 particles traverse 
the tile elements. Two different event simulators, HIJING [13] and Fritiof [14], 
were used to check the sensitivity of the analysis procedure to the assumed pri- 
mary particle distributions and the corresponding momentum distributions. 
The deduced particle multiplicities using the different codes agree to better 



3 



than 3%. This uncertainty is folded into the quoted overall systematic uncer- 
tainties. The particle yields are corrected for background contributions using 
the GEANT simulations of the array response. These corrections are position 
dependent and range from 20% to 40% for the TMA and from 6% to 25% for 
the SiMA, increasing with increasing \r)\. 

The BBC Arrays consist of two sets of Cherenkov UV transmitting plastic 
radiators coupled to photomultiplier tubes. They are positioned around the 
beam pipe on either side of the nominal interaction point at a distance of 
2.15 m. One array consists of 8 "large" detectors with 51 mm diameter pho- 
totubes and 36 "small" detectors with 19 mm diameter phototubes, arranged 
symmetrically around the beam pipe. The other array is asymmetric to allow 
the movement of the forward spectrometer to small angles and consists of 5 
large detectors and 30 small detectors. The BBC elements have an intrinsic 
time resolution of 65 ps allowing for the determination of the position of the 
interaction point with a precision of as 1.6 cm. Charged particle multiplicities 
in the pseudorapidity range 2.1 < \rj\ < 4.7 are deduced from the number of 
particles hitting each tube, as found by dividing the measured ADC signal by 
that corresponding to a single particle hitting the detector. Corrections for 
background events were based on GEANT simulations, with the correction 
factors ranging from 37% to 50% of the measured yield, depending on vertex 
location. 

The two ZDCs, positioned on either side of the nominal interaction point at 
a distance of 18 m, were used to establish a minimum bias trigger for the 
experiment. The ZDCs, which consist of alternating layers of tungsten and 
fiber-ribbon plastic scintillator, measure the energy deposited by spectator 
neutrons that are emitted at small angles with respect to the beam direction, 
with \6\ < 2 mrad [10]. The time difference between the two ZDC signals can 
be used to determine the interaction vertex with a precision of ~ 3.6 cm. 

In addition to the global detector measurements, charged particle multiplicities 
were also obtained at 40°, 60°, and 90°, corresponding to 77=1.1, 0.55, and 
0.0, respectively, using the TPM1 counter. This detector, which is operated 
in a field-free region, has an active volume of (w,h,d) = (37.5 cm, 21.0 cm, 
36.0 cm) with its volume center located 94.5 cm from the nominal vertex. For 
TPM1, particle multiplicities are deduced by counting reconstructed tracks 
arising from collision vertices within 15 cm of the nominal vertex position. The 
acceptance in 77 for this counter, accounting for both the geometric acceptance 
and the range of vertex locations, varies from Ar/=0.5 at 90° to A??=0.60 at 
40°. Through the track reconstruction, the interaction vertex can be located 
to better than 0.5 cm. A correction for the background of secondary particles 
was found by projecting the vertical position of tracks in TPM1 back to the 
beamline. The resulting spectrum was fitted with the sum of three Gaussian 
distributions; one for the peak of primary particles and summed background 
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contributions from weak decays or secondary interactions, and delta electrons 
from the beampipe, respectively. The ratio of the integral of the signal peak 
to the two background distributions was taken as the fraction of primary 
particles. Typical background-to-total ratios of ~5% were found at all angles. 

The dominant systematic uncertainties for the TPM1 measurements arise 
from the background subtraction and a tracking efficiency correction. The 
efficiency correction is obtained by inserting simulated single tracks into ac- 
tual data events and analyzing the reconstruction probability using the same 
track reconstruction code employed for the experimental results. The track- 
ing efficiency in the TPM1 counter ranges from almost unity for low density, 
peripheral scattering events, to ~0.9 for central events at 90° and ~0.85 for 
central events at 40°. In contrast, the dominant systematic uncertainty near 
mid-rapidity for the SiMA and TMA arrays reflect the single-particle response 
calibrations obtained for the constituent detectors of these arrays. At larger 
pseudorapidity values, the uncertainty in the background subtraction also be- 
comes significant. In addition, the relatively thick scintillator elements of the 
TMA array are expected to be somewhat sensitive to beam related background 
that does not come from the collision vertex. These effects lead to the some- 
what larger uncertainty quoted for the TMA as compared to the SiMA. The 
BBC systematic uncertainty is dominated by the background correction. Over- 
all, the systematic uncertainties for the different detector systems are: SiMA, 
±8% for 77<1.5, increasing to ±10% for rj >2.5; TMA, ±12% for rj <1.5 in- 
creasing to ±15% for rj >2.0; BBC, ±10%; TPM1, ±9%/-7% for the most 
central events and ±6% for the most peripheral. 

The centrality selection for the experiment was obtained by developing a 
minimum-biased multiplicity distribution using the MA, reflecting all events 
for which there is a nuclear interaction, and assuming that a cut on the total 
multiplicity translates to a cut on collision centrality. Fig. 1 shows the multi- 
plicity distribution of charged particles established by the MA, normalized to 
the maximum observed multiplicity. The independent multiplicity measure- 
ments of the SiMA and TMA detectors are summed for this figure, after the 
SiMA multiplicity was rescaled to account for the difference in the geometric 
coverage of the two arrays. Where possible, the location of the interaction 
vertex was determined using TPM1. Otherwise, the vertex determined from 
the BBC or ZDC data was adopted. The charged particle densities reported in 
this paper are obtained almost exclusively using the TPM1 and BBC vertices 
and are relatively insensitive to the specific choice. The multiplicity distribu- 
tion requires that a coincidence occurs between the two ZDC detectors, an 
interaction vertex be located within 30 cm of the array center (as obtained by 
one of the three possible vertex measurements), and that a TMA multiplic- 
ity of at least four be observed. The TMA multiplicity cut eliminates events 
corresponding to beam-gas interactions and for cases where tails of pedestals 
for individual detector elements add to simulate real collision events. Based 
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on HIJING simulations, where the shape of the multiplicity distribution for 
events with fewer than 200 tile hits is used to estimate the missing yield for 
events with fewer than four tile hits, it is estimated that the TMA array ob- 
serves cleanly 95% of the total nuclear cross sections. This value is used to 
correct for the missing cross section where the tile multiplicity is less than four 
when determining the final multiplicities. 

Centrality selection in the BBC analysis was done using the multiplicity dis- 
tribution obtained by summing the hits in the BBC arrays. This allows for the 
use of interaction vertices within 120 cm of the MA center, well outside the 
range for which the MA multiplicities could be reliably determined. The insert 
in Fig. 1 shows the strong correlation of the BBC and MA multiplicities, nor- 
malized to their maximum values. In the pseudorapidity range of 3.0< r] <4.2, 
where it was possible to analyze the BBC data using both centrality selections, 
the two analyses give identical results to within 1%. Still, the lower segmen- 
tation of the BBC arrays leads to a less precise multiplicity measurements as 
compared to the MA. 

Figure 2 presents the measured dN ch /dr) distributions for a number of central- 
ity cuts. The points associated with positive pseudorapidity values correspond 
to detectors located on the forward spectrometer side of the BRAHMS appa- 
ratus. Within the quoted systematic uncertainties, good agreement is found 
among the measurements using different detector systems. There is, however, 
an offset observed for the TMA data which may reflect non-collision related 
background events. The other detectors are relatively insensitive to such a 
background mechanism. 

Particle densities at selected pseudorapidities are tabulated in Table 1. This 
table also presents the integrated charge-particle multiplicities in the range 
—4.7 < rj < 4.7. For central events at mid-rapidity, the observed particle 
densities are about 1.8 times greater than observed for central events in Pb+Pb 
scattering at y/s NN =17.2 GeV [4]. 

The PHOBOS [6] and PHENIX [8] experiments have recently reported values 
for the charged particle densities at mid-rapidity corresponding to the 6% 
most central events. PHOBOS obtains dN c h/di]\ v=Q = 555±l2(stat) ±35(sy st) 
while PHENIX finds dN ch / 'c% =0 = 609 ± l(stat) ± 37 '(syst). For the same 
centrality range, our result is dN c h/ d7]\ v= o = 549 ± l(stat) ± 35(syst). Within 
the systematic uncertainties the three experiments are in agreement. 

Fig. 3 shows the charged-particle multiplicity distributions scaled by half of 
the average number of participant nucleons for the 0-5% and 30-40% central- 
ity cuts. The average number of participants (< N part >) corresponding to 
the different centrality cuts was obtained using the HIJING code [13], where 
the collision geometry is based on the Glauber model [15]. The uncertainty in 
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< N par t > largely reflects the uncertainty in the multiplicity measurement and 
is estimated as <0.5% for 5% centrality, increasing to 4% at 50% centrality. 
For Fig. 3 (and Fig. 4), the dN ch /dr] distributions have been symmetrized and 
uncertainty weighted averages have been developed for the different global 
detector subsystems. It is found that near mid-rapidity, the charged-particle 
yield per participant pair is significantly greater for more central events as 
compared to peripheral events. However, at larger pseudorapidities, this dif- 
ference disappears. The mid-rapidity behavior is consistent with that observed 
by the PHENIX [8] and PHOBOS [7] collaborations and attributed [8] to an 
increased contribution of harder scattering processes for more central colli- 
sions. 

The insert in Fig. 3 shows the Au+Au results for large pseudorapidities in 
terms of the pseudorapidity shifted by the beam rapidity, thus allowing for 
a comparison with SPS Pb+Pb data [4] within the "limiting fragmentation" 
picture [16]. This hypothesis states that at high energies the number of par- 
ticles produced by the "wounded projectile" nucleons should be independent 
of the details of the "target", the "projectile", and the beam energy. When 
shifted by the beam rapidity, effectively translating to the beam's reference 
frame (assuming similar values for the pseudorapidy and rapidity variables), 
the normalized yields near r}'{— rj — ybeam) = should be system and energy 
independent. A "limiting fragmentation" behavior has been observed in pp, 
pp, p-emulsion, and 7r-emulsion data (see references 34-37 of ref. [4]). It is 
shown in ref. [4] that the hypothesis also works in comparing Pb+Pb data at 
v /s^7=17.2 GeV with 200 GeV/nucleon O+AgBr and S+AgBr results [17]. 
The closed triangles in the Fig. 3 insert show the Pb+Pb results. Inspection 
of this figure, where the current results at v /J/v r N r =130 GeV follow the same 
general curve as found for the lower energy data, suggests that the "limiting 
fragmentation" picture is successful over a very large range of energies ex- 
tending from SPS to RHIC ( v /5^=17.2 GeV to 130 GeV). The absence of a 
centrality dependence for the scaled yields near beam rapidity, as seen in Fig. 
3, is also consistent with such a picture. 

In Fig. 4 we compare HIJING [13], UrQMD [18], and AMPT model [19-21] 
calculations to the data for different centrality ranges. HIJING is a model 
that places special emphasis on perturbative QCD processes leading to mul- 
tiple minijet production from parton scatterings. The model, which does not 
include final state rescattering except for schematic jet quenching, is found to 
do a good job in describing the data at mid-rapidity, but underestimates the 
widths of the multiplicity distributions at all centralities. The AMPT model 
uses HIJING to generate the initial phase space of partons, but then extends 
the calculations to model the parton-parton collisions and the final hadronic 
interactions. By including these final state interactions, the widths of the mul- 
tiplicity distributions are found to agree well with experiment, suggesting the 
importance of these final state interactions. It is also interesting to note that 
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the LUCIFER code [22] gives similar results to AMPT but within a purely 
hadronic formulation. The UrQMD transport model accounts for hadronic in- 
teractions at low and intermediate energies in terms of interactions between 
known hadrons and their resonances and at higher energies in terms of exci- 
tations of color strings and their subsequent fragmentation. For the present 
RHIC data, the shape of the pseudorapidity distributions are well reproduced, 
however with the particle production being way too large. The current results 
appear to rule out the UrQMD model as a description of the RHIC data. 

The BRAHMS experiment has measured pseudorapidity densities of charged 
particles from Au+Au collisions at v /s/v r Jv=130 GeV over a large range of 
pseudorapidity as a function of collision centrality using several independent 
detector systems and methods. We observe an enhancement of particle pro- 
duction for central collisions at mid-rapidity, suggesting a breakdown in the 
simple scaling by the number of participant nucleons seen at larger pseudo- 
rapidities. Comparisons of the Au+Au fragmentation region data with SPS 
energy Pb+Pb results suggest that a "limiting fragmentation" behavior is 
seen for nucleus- nucleus collisions and is, in fact, already reached at the lower 
energy. With the inclusion of final hadronic interactions, parton scattering 
models are found to give a good description of the observed behaviors. 
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Table 1 

Charged particle densities in dN c h/dri as a function of centrality and pseudorapid- 
ity. Total uncertainties, dominated by the systematics, are indicated. The average 
number of participants < N part > is given for each centrality class based on HIJING 
model calculations. The last column gives the integral charged particle multiplicity 
within the pseudorapidity range —4.7 < 77 < 4.7. 



Centrality 


< N, part > 


77 = 


7] = 1.5 


r] = 3.0 


T] = 4.5 


N ch 


0-5% 


352 


553±36 


554±37 


372±37 


107±15 


3860±300 


5-10% 


299 


447±29 


454±31 


312±36 


94±13 


3180±250 


10-20% 


235 


345±23 


348±25 


243±27 


79±10 


2470±190 


20-30% 


165 


237±16 


239±16 


172±18 


59±8 


1720±130 


30-40% 


114 


156±11 


159±11 


117±13 


43±6 


1160±90 


40-50% 


75 


98±7 


104±7 


77±9 


30±4 


750±60 
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Fig. 1. Normalized MA array multiplicity distribution, as discussed in the text. 
The scatter plot insert shows the correlation of the normalized BBC and MA mul- 
tiplicities. A schematic drawing of the global detectors used in the measurement, 
excluding the ZDCs, is also shown. 



11 



600 
400 

— 200 

£°100- 
^ 80 

60 
40- 

20 



a/a* gg8@S9@^§tS^@g @n a_\A 



6 



a tr 

// 

/ 



*V a' 
a/ A 

y* , 



□□□□□□□□□□ 

oooooqog^oc^oooooo 



A" 



A A 

A 



/ 



o SiMA 
□ TMA 

a A BBC 
*TPM1 

_j i i l_ 



A 

\ A 

A 

\ 

_l ^_ 



V A 



-4 -2 2 4 

n 

Fig. 2. Distributions of dN c h/dr] for centrality ranges of, top to bottom, 0-5%, 
5-10%, 10-20%, 20-30%, 30-40%, and 40-50%. Statistical uncertainties are shown 
where larger than the symbol size. The connecting lines and alternating open and 
closed symbols for the BBC data are to help distinguish points associated with 
different centrality ranges. 
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Fig. 3. Charged particle densities normalized to the number of participant pairs. 
Total (statistical + systematic) uncertainties are shown, not including the uncer- 
tainty in < Npart > (see text). The insert compares the 0-5% central (open circles) 
and 30-40% central (open squares) Au+Au results to the 9.4% central Pb+Pb data 
(closed triangles) of ref. [4]. For this comparison of the "projectile" regions for the 
two reactions, the data are plotted in terms of the pseudorapidity shifted by the 
beam rapidity, as discussed in the text. For the insert, the particle densities are 
normalized to the number of projectile participants, which is equal to the number 
of participant pairs for a symmetric reaction. 
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Fig. 4. Distribution of dN c h/dr] for the indicated centrality ranges. Total uncertain- 
ties are indicated. HIJING, AMPT, and UrQMD model calculations are shown. 
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